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ABSTRACT 

Using high-resolution (R ~ 22,000) near-infrared (1.51 - 1.75 /jm) spectra from Keck Observatory, we 
measure the kinematic masses of two super star clusters (SSCs) in M82. Cross-correlation of the spectra 
with template spectra of cool evolved stars gives stellar velocity dispersions of a r — 15.9 ± 0.8 km s _1 
for J0955505+694045 ('MGG-9') and oy = 11.4 ± 0.8 km s" 1 for J0955502+694045 ('MGG-11'). The 
cluster spectra are dominated by the light of red supergiants, and correlate most closely with template 
supergiants of spectral types M0 and M4.5. King model fits to the observed profiles of the clusters in 
archival HST/NICMOS images give half-light radii of r hp = 2.6 ±0.4 pc for MGG-9 and r hp = 1.2 ±0.17 
pc for MGG-11. Applying the virial theorem, we determine masses of 1.5±0.3 x 10 6 M Q for MGG-9 and 
3.5 ± 0.7 x 10 5 M Q for MGG-11 (where the quoted errors include oy, rn p and the distance). Population 
synthesis modelling suggests that MGG-9 is consistent with a standard initial mass function, whereas 
MGG-11 appears to be deficient in low-mass stars relative to a standard IMF. There is, however, evidence 
of mass segregation in the clusters, in which case the virial mass estimates would represent lower limits. 

Subject headings: galaxies: individual (M82) — galaxies: starburst — galaxies: star clusters — 
galaxies: nuclei — infrared: galaxies 



1. INTRODUCTION 

Starburst galaxies are the site of about 25% of the high- 
mass (> 8M Q ) star formation in the local Universe (Heck- 
man 1998). As high-mass stars form predominantly in 
dense clusters and OB associations (Miller & Scalo 1978), 
massive stellar clusters in nearby starburst galaxies pro- 
vide a laboratory for studying intense star formation and 
related feedback processes. Star formation in starbursts is 
resolved into massive, compact super star clusters (SSCs) 
which account for perhaps 20% of the new stars formed in 
the burst (Meurer et al. 1995; Zepf et al. 1999). The combi- 
nation of high- resolution HST imaging and high- resolution 
near-infrared and optical spectroscopy has found SSC sizes 
(~ 1-5 pc) and masses (10 5 -10 6 M ) comparable with 
those of globular clusters. Estimates of cluster ages have 
shown the SSCs to be much older than the crossing time, 
indicating that they are gravitationally bound. In addi- 
tion, radial profiles of the clusters are best fit with the 
King models (e.g., Mengel et al. 2002) derived for globu- 
lar clusters (King 1962). These observations suggest that 
SSCs are young globular clusters, forming in the present- 
day universe. 

While SSCs have been identified in a number of galaxies, 
kinematic masses have been measured for relatively few. 
Among these are dwarf galaxies NGC 1569 (Ho & Fil- 
ippenko 1996a; Gilbert & Graham 2001) and NGC 1705 
(Ho & Filippenko 1996b), and the NGC 4038/4039 merger 
(Mengel et al. 2002). While much of this work has relied 
upon optical spectroscopy, there are distinct advantages 
in the near-infrared. SSCs are often situated in dusty re- 
gions where extinction is high; the most heavily-reddened 
clusters are not seen in optical images. In addition, the 



spectra of SSCs older than ~ 7 Myr are dominated by the 
light of cool supergiant stars. The near-IR spectra of these 
evolved stars feature a multitude of atomic and molecu- 
lar absorption lines which are suitable for measurement of 
the stellar velocity dispersion within a cluster via spectral 
cross-correlation. Combined with high-resolution images, 
this technique allows determination of cluster masses. 

Kinematic masses for SSCs are relevant to the debate 
over the form of the initial mass function (IMF) in star- 
burst galaxies (e.g., Eisenhauer 2001; Gilmore 2001). No- 
tably, Rieke et al. (1993) use population synthesis models 
for the nuclear starburst in M82 to conclude that the IMF 
is deficient in low- mass stars (< 3M Q ) relative to the IMF 
of the solar neighborhood. Sternberg (1998) finds that 
SSC NGC 1705-1 appears to be deficient in low-mass stars, 
whereas a large fraction of the stellar mass in the SSC NGC 
1569A is contained in low-mass stars. Measurement of the 
kinematic mass of an SSC provides the only means of char- 
acterizing the low-mass stellar population of the cluster. 
While the light of these stars is overwhelmed by the far 
more luminous supergiants, the low-mass stars comprise a 
substantial fraction of the cluster's mass for a normal IMF. 
Quantifying the low-mass stellar content is critical to con- 
straining the IMF as well as assessing whether a cluster 
can remain bound as a globular cluster. 

M82 (NGC 3034) provides an interesting case. As one 
of the nearest starburst galaxies (3.6 Mpc) it has been 
extensively observed from X-rays to radio wavelengths. 
Its proximity results in an obvious resolution advantage: 
star-forming regions can be studied on small spatial scales 
(1" = 17.5 pc) and individual SSCs are resolved by HST. 
A NICMOS F160W mosaic (Figure 1) of the nuclear star- 
burst region shows many luminous SSCs within ~ 300 pc 
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of the nucleus. The nuclear starburst is thought to be 
active in the sense that the typical age for the starburst 
clusters is <~ 10 7 years (Satyapal et al. 1997). Evolution- 
ary synthesis models by Forster Schreiber (1998) suggest 
the nuclear starburst (i.e., the central 450 pc) consists of 
two distinct events with ages of about 5 Myr and 10 Myr. 
O'Connell et al. (1995) imaged M82 in the V- and /-bands 
with the high-resolution Planetary Camera aboard HST, 
identifying over 100 SSCs within a few hundred parsecs 
of the nucleus, de Grijs et al. (2001) imaged a region in 
the disk of M82, 1 kpc from the nuclear starburst, with 
WFPC2 and NICMOS. They identified 113 SSC candi- 
dates which they believe were part of a starburst episode 
<~ 600 Myr ago (a "fossil starburst"). Smith & Gallagher 
(2001) derived the kinematic mass of SSC 'M82-F' us- 
ing high-resolution optical (6010-8132A) spectroscopy and 
WFPC2 (F439W) images. M82-F (hereinafter MGG-F) 
lies ~ 500 pc west of the nucleus of M82; Smith & Gal- 
lagher estimate an age of 60 ± 20 Myr, intermediate be- 
tween the ongoing nuclear burst and the fossil burst farther 
out in the disk. 

In this paper, we investigate the kinematic mass of 
two clusters (J0955502+694045, hereinafter MGG-11, and 
J0955505+694045, hereinafter MGG-9) on the western 
edge of the active nuclear starburst. These particular 
clusters were chosen for this pilot program based on their 
brightness (run ~ 13) and the ability to place both on a 
single slit position. 

Frcedman et al. (1994) measured the distance to neigh- 
boring galaxy M81 using Cepheids, determining a value of 
3.6 ± 0.3 Mpc. Observations of HI (e.g., Applcton et al. 
1981) indicate a physical link between M81 and M82. As 
the projected separation between the galaxies is only <~ 40 
kpc, well within the uncertainty on the distance, we adopt 
this value for M82. The corresponding distance modulus 
ism-M = 27.8. 

2. OBSERVATIONS 

We observed super star clusters MGG-9 and MGG-11 
with the 10-m Keck II telescope on 2002 February 23, 
using the facility near-infrared echelle spectrometer NIR- 
SPEC (McLean et al. 1998). We obtained high-resolution 
(R ~ 22,000), cross-dispersed spectra in the wavelength 
range 1.51-1.75 fim using the NIRSPEC-5 order-sorting 
filter. The data fall in seven echelle orders, ranging from 44 
through 50. The clusters were imaged in four nods along 
the 0."432 x 24" slit, with successive nods separated by 5" 
(spatial resolution of the spectra is 0."5). Each integration 
was 600 seconds, providing a total time of 2400 seconds on 
the clusters. The spectra were dark subtracted, flat-fielded 
and corrected for cosmic rays and bad pixels. The curved 
echelle orders were then rectified onto an orthogonal slit- 
position/ wavelength grid based on a wavelength solution 
from sky (OH) emission lines. Each pixel in the grid has a 
width of S\ = 0.19 A. We sky-subtracted by fitting third- 
order polynomials to the 2D spectra column-by-column. 
The individual extracted spectra had average signal-to- 
noise ratios of - 47 for MGG-9 and - 42 for MGG-11. 

The NIRSPEC echelle turret is jostled whenever the 
cryogenic image rotator is slewed. When large slews, 
with amplitudes comparable to the permitted range of mo- 
tion, are executed, the echellegram shifts by several pixels. 



However, when the image rotator is making the small, slow 
motions necessary to remove sidereal rotation, the wave- 
length solution is stable to better than a few hundredths 
of a pixel. This is true even when an object is crossing 
the meridian and the sidereal rotation is fastest. Thus the 
velocity broadening reported here is intrinsic to the source 
and not an instrumental effect. 

The cluster spectra were extracted using Gaussian 
weighting functions matched to the wavelength-integrated 
profiles. To correct for atmospheric absorption, we ob- 
served a B8V star, HD 74604. Both the clusters and the 
calibration star were observed at an airmass of ~ 1.6. To 
account for photospheric absorption features (particularly 
Brackett and helium lines) and continuum slope, the cali- 
bration star spectrum was divided by a spline function fit. 
The resulting atmospheric absorption spectrum was then 
divided into the cluster spectra. 

We also observed several late-type stars of luminosity 
class I, II or III with NIRSPEC for use as spectroscopic 
templates in cross-correlation analysis (see Section 3). Ta- 
ble 1 summarizes these observations. Template stars were 
observed with the same slit, filter, and cross-disperser and 
echelle grating positions as the clusters. Reduction, recti- 
fication and extraction steps were identical to those de- 
scribed above for the clusters, with the exception that 
sky subtraction for the high signal-to-noise stars was per- 
formed by pairwisc subtraction of successive nods along 
the slit. 

For photometry and determination of structural param- 
eters of super star clusters in M82, we obtained images 
from the HST archive (HST Proposal 7218, Marcia Rieke, 
PI). The images were taken 1998 April 12 with the NIC2 
camera on NICMOS and cover the nuclear region with four 
adjacent, overlapping fields along the major axis. We used 
the images taken with the broadband F160W and F222M 
filters. NIC2 has a plate scale of 75 mas per pixel, and 
is therefore critically sampled at 1.75/im. Each field was 
observed in a spiral dither of four images, resulting in a 
field of view of <~ 19. "5 x 19. "5 for each pointing center. 
We relied upon the NICMOS data reduction pipeline (cal- 
nicb) for combination of the dithered images. The dither 
pattern used a point-spacing of 0."215 (2.87 pixels). At 
the distance of M82, one pixel corresponds to 1.3 pc. A 
mosaic of the four fields covers an area of approximately 
29" x 72". Total exposure time for each of the four fields 
was 192 seconds in F160W and 448 seconds in F222M 
(except the western-most field, which was 576 seconds). 
We relied on the standard pipeline reduction procedures 
(calnica) for dark-subtraction, flat-fielding and cosmic ray 
rejection, and did not perform any additional calibration 
steps ourselves. 

We compared the astrometry of identifiably common 
clusters between the NICMOS images and the positions 
listed in Kronberg et al. (1972), and found discrepancies 
of ~ 2". All astrometry quoted in this paper has been 
adjusted to match Kronberg et al., and is correspondingly 
accurate to ~ 0."6. 

3. CLUSTER VELOCITY DISPERSION 

In the near-infrared, SSC spectra are dominated by the 
light of cool evolved stars, the spectra of which contain 
a large number of molecular and atomic absorption fea- 
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tures (Gilbert & Graham 2001). We have assembled high- 
resolution NIRSPEC spectra of a number of evolved (lu- 
minosity class I, II and III) stars, ranging in spectral type 
from G2 through M4.5. Cross-correlation of these tem- 
plate spectra with the spectrum of an individual super 
star cluster allows identification of the dominant spectral 
type of stars in the SSC as well as determination of the 
line-of-sight velocity dispersion. For each cluster in the 
present work, we have 28 measurements to cross-correlate 
with the template star spectra: four nods in each of seven 
echelle orders. We have chosen not to sum the spectra for 
a given echelle order to avoid the interpolation required 
in registering the separate measurements to a consistent 
wavelength grid and the resulting smearing of the spectral 
resolution. 

Certain atmospheric OH emission lines were incom- 
pletely removed in the sky subtraction process. The resid- 
uals of these bright lines, generally close doublets, must be 
removed prior to cross-correlation to avoid introduction of 
systematic bias. We have replaced the pixels affected by 
residual sky emission with the median value of the five 
pixels on either side of the contaminated range. The frac- 
tion of pixels replaced in a given echelle order typically 
amounts to a few percent. Tests on an OH emission-free 
echelle order indicate that the cross-correlation result is 
unaffected within the stated uncertainities, with the fol- 
lowing exception. At the systemic velocity of M82, the 
bright telluric OH 4 — 2 P2(3.5) emission line at 1.6194/im 
overlies the CO (v — 3 — 6) bandhead in echelle order 47 
(1.609/um < A < 1.636/xm). The CO bandhead feature is 
so prominent in the spectrum that residual noise from sky 
subtraction introduces significant systematic error for this 
echelle order. The results of the cross-correlation analysis 
for echelle order 47 are therefore unreliable, and have not 
been used in this analysis. 

The remaining 24 spectra for each cluster are cross- 
correlated with the corresponding echelle orders of the 
template stars, and the template star spectra are auto- 
correlated. Our cross-correlation is based on the method 
of Tonry & Davis (1979). We fit a fifth-order Chebychcv 
polynomial to the spectrum to subtract the continuum, 
and apodize with a cosine bell to remove mismatch at 
the edges. The spectrum is then transformed into Fourier 
space, and multiplied by fourth-order Butterworth filters 
to remove low-frequency continuum variations and high- 
frequency noise. The filter frequencies are selected to 
minimize any systematic bias introduced by the presence 
of noise. By the convolution theorem, multiplication (in 
Fourier space) of the continuum-subtracted, apodized, fil- 
tered spectrum and its complex conjugate, followed by in- 
verse Fourier transformation, gives the cross-correlation 
function. Invoking the central limit theorem, we fit the 
auto- and cross-correlation peaks with Gaussians to de- 
termine their respective widths. The excess width of the 
cross-correlation relative to the auto-correlation of a tem- 
plate star is a measure of the stellar velocity dispersion 
within the cluster. Cross-correlation of a pair of template 
stars should reflect only differences in rotation (vsini), 
macroturbulence, and noise. The advantage of using the 
auto-correlation is that it fully accounts for the intrinsic 
width of absorption features when the correct template 
spectral type is used. Only the dispersion of photospheric 



macroturbulence within that spectral type and luminos- 
ity class affects our analysis, and no first-order correction 
to the measured velocity dispersion is necessary. The ob- 
served macroturbulence for K supergiants is 11 km s _1 
with a peak-to-peak range of ~ 2 km s _1 (Gray & Toner 
1987). Our analysis assumes that our template stars are 
representative of their spectral type and luminosity class. 

One potential difficulty is metallicity differences between 
Galactic supergiants used as templates and the supergiants 
producing the cluster light. McLeod et al. (1993, and refs.) 
cited evidence from emission-line studies of various ele- 
ments and concluded that the present-day ISM in M82 has 
solar or slightly higher metallicity. Forster Schreiber et al. 
(2001) determined that near-IR stellar absorption features 
observed in the starburst core are consistent with the light 
from solar- metallicity red supergiants (RSG). Metallicity 
factors arc unlikely to significantly bias our measured ve- 
locity dispersions. 

The results of the cross-correlation with the template 
stars are presented in Table 2. The peak amplitude of 
the cross-correlation function (CCF) measures the simi- 
larity of the cluster to the template spectral type. The 
CCF peak increases steadily with spectral type from 0.3 
for the G2I template to 0.8 for the M0I and M4.5I tem- 
plates. We find that the spectra of clusters MGG-9 and 
MGG-11 closely match the spectrum of the MOIab (e.g., 
Figures 2 and 3) and M4.5Iab stars. The average of the 
CCF over the 24 spectra for each cluster is shown in Fig- 
ure 4. Unfortunately, our spectral atlas does not presently 
include any template supergiant spectra between M0 and 
M4.5. Intermediate spectral types, e.g. M2, should match 
the clusters at least as well as M0 and M4.5. 

Asymmetries in the CCF indicate kinematic substruc- 
ture, such as the presence of an additional cluster. NGC 
1569A, for example, displays an asymmetric CCF which 
is the result of a kinematically distinct cluster which is 
projected onto the same slit position (Gilbert & Graham 
2002; Ho & Filippenko 1996a). The CCFs for these M82 
SSCs are quite symmetrical, and are well fit by Gaussians. 
The clusters show evidence of only a single kinematic com- 
ponent. 

Based on an average of the results from the best-match 
templates, we measure velocity dispersions of oy = 15.9 ± 
0.8 km s- 1 for MGG-9 and a r = 11.4 ± 0.8 km s" 1 for 
MGG-11. (The uncertainties quoted here are the stan- 
dard errors — the standard deviation of the independent 
measurements divided by y/N.) 

Cross-correlation with a mismatched template spectrum 
can introduce systematic bias to the velocity dispersion 
determination. We note a roughly linear trend between 
the measured velocity dispersion and the peak value of 
the CCF for a given template. In the case of MGG-11, 
extrapolation of the trend to a hypothetical "perfect" cor- 
relation of 1.0 gives a value of ay which is consistent with 
the value quoted above. In the case of MGG-9, the strong 
decreasing trend of oy with the CCF peak value suggests 
that the quoted value may be slightly inflated due to sys- 
tematic bias from template mismatch. Comparison to a 
correlation of 1.0, however, is not the correct metric; the 
velocity dispersion of the cluster stars reduces the CCF 
peak even for a perfect spectral match. As a test, we 
convolved the MOIab spectrum with a Gaussian velocity 
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distribution, and then cross-correlated the result with the 
template MOIab spectrum. The CCF peaked at 0.90 for 
a velocity dispersion of 16.2 km s _1 ; increasing the veloc- 
ity dispersion decreases the CCF peak (note the higher 
CCF peak values for MGG-11, which has a lower velocity 
dispersion than MGG-9). We adopt the quoted ay values 
with the caveat that there may be a small amount (<~ 1 km 
s _1 ) of systematic bias present due to template mismatch. 

Cross-correlation analysis also measures the heliocentric 
radial velocities of the clusters. We find that MGG-9 has 
v r = 113.8 ±0.4 km s^ 1 and MGG-11 has v r = 92.5 ±0.5 
km s _1 . These values are based on cross-correlation with 
HD 43335, which has a radial velocity of 38.9 km s -1 (per 
SIMBAD 3 ). The spectrum of HD 43335 is was taken imme- 
diately before the cluster spectra, and we find no grating 
shifts which would systematically bias the radial velocities. 

Millimeter observations of CO emission by Shen & Lo 
(1995) provide context for the radial velocities. The clus- 
ters lie in a "saddle" feature in the velocity-integrated CO 
emission map, roughly 10" west of the kinematic center. 
Plotted on Shen & Lo's major axis p-v cut, the clusters 
lie on the western shoulder of a prominent second spatial 
component, separate from the primary velocity gradient. 
The radial velocities of the clusters measured in the near- 
IR are consistent with the velocity of CO at their locations, 
suggesting a potential kinematic link between the clusters 
and their natal molecular gas. Dcprojecting the galaxy 
and applying the Shen & Lo CO velocity gradient to a 
simple cylindrical solid-body rotation model gives an or- 
bital period of ~ 17 Myr for the clusters. They have not 
had sufficient time to complete an orbit of the galaxy, and 
may yet be near their birthplaces. Future high-resolution 
studies will shed more light on the kinematic connection 
between dense molecular gas and young super star clusters 
in the starburst. 

4. CLUSTER INTEGRATED LUMINOSITY 

Figure 1 shows the NICMOS F160W mosaic of the nu- 
clear starburst in M82, with the clusters treated in this pa- 
per labelled with our secondary designations. The region 
is notoriously complex, showing structure on virtually any 
spatial scale resolvable. In the NICMOS images, many of 
the candidate super star clusters appear against a mottled 
background of relatively strong diffuse emission. In order 
to perform photometry and measure size parameters, we 
must isolate the cluster light from the background. We 
placed a 0."15-wide annulus at an inner radius of 0."5 (9 
pc) from the centroid of each cluster, and fit a minimum- 
least-squares plane to the points within the annulus. Sub- 
traction of this plane removes the background emission, 
allowing us to fit the cluster light. 

Calculation of the virial mass of a cluster (Section 5.1) 
requires determination of the half-light radius. Special 
care is required, in that many of the clusters are only 
marginally resolved. The effect of the instrumental point 
spread function must therefore be taken into account. For 
each cluster, we generated a model NIC2 PSF using the 
Tiny Tim program (Krist 1995). The PSF models were 
matched to the centroid location of each cluster on the de- 
tector, oversampled 2x and rebinned to match the pixel 
size of NIC2. To model the observed cluster images, we 

3 http://simbad.u-strasbg.fr/ 



tested various functional forms of the radial brightness pro- 
file. We found good fits to the data for the empirical King 
model (King 1962) of the form: 



\[l + (r/r c ) 2 ]5 [l + (r t /r c ) 2 ]5 J 

where < r < r t . The King model has three free param- 
eters: the core radius r c , the central concentration r t /r c 
(where r t is the tidal radius), and a scaling factor k. The 
King model is convolved with the PSF model and com- 
pared to the data. The best-fit model is determined us- 
ing a Levenberg-Marquardt minimization, finding the least 
squares fit to the data through iterative search of param- 
eter space. The resulting models give very good fits, as 
illustrated in Figure 5. Once the King model parameters 
have been fit for a given cluster, the projected half-light 
radius (r^ p ) is the numerical solution of 

J^ hp f(r)2irrdr _ 1 
So* f(r)2irrdr 2 

Likewise, integration of the fitted King model over the 
interval < r < r t , in conjunction with the NIC2 photo- 
metric keyword PHOTFNU, provides the total apparent 
magnitude of the cluster. Table 3 lists the astrometry, pho- 
tometry and derived half-light radii for the 20 NICMOS- 
resolved clusters labelled in Figure 1. In particular, we find 
projected half-light radii (r^) of 146 ± 16 mas for MGG-9 
and 66 ± 7 mas for MGG-11. At the adopted distance of 
3.6 ± 0. Mpc, the projected half-light radii are 2.6 ± 0.4 pc 
for MGG-9 and 1.2 ± 0.17 pc for MGG-11. 

To determine the uncertainties associated with the pho- 
tometry and half-light radii, we performed a Monte Carlo 
simulation wherein we added fake star clusters to the NIC- 
MOS images and used our fitting algorithm to recover 
them. We tested fake clusters with a range of luminosities 
and King model parameters, placed against a representa- 
tive range of background conditions. While the algorithm 
could not reliably recover the values of r c and central con- 
centration (r c /r t ), it must be noted that these parameters 
are highly covariant; a particular halflight radius may be 
represented by a family of r c and rt values. The fitting al- 
gorithm identified the projected halflight radii to within an 
uncertainty of 12%, independent of the input values of r c 
and r t . Photometric uncertainities scaled with cluster lu- 
minosity, ranging from 0.03 to 0.4 mag in low background 
conditions. In regions of high and/or variable background 
emission, the algorithm produced larger scatter around the 
correct luminosity. The results of this analysis are reflected 
in the photometric uncertainties in Table 3. 

We have chosen to use the F160W images for determina- 
tion of the half-light radii rather than the F222M images 
for two reasons: (1) we are measuring the velocity disper- 
sion of our test particle stars at 1.51 - 1.75 fim, and it is 
consistent to determine structural parameters of the test 
particles by using the same light; and (2) warm circumstel- 
lar dust may result in excess continuum emission longward 
of 2 /im, complicating cross-correlation analysis. 

MGG-F has been closely examined by other au- 
thors, and provides a comparison for our measurements. 
O'Conncll ct al. (1995) imaged MGG-F in the V-band 
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(F555W) with WFPC. They deconvolve the images using a 
Tiny Tim PSF model, and find a mean FWHM of 160 mas 
(2.8 pc). The authors note that the deconvolution is "not 
highly reliable" because the cluster falls near the edge of 
the detector. Smith & Gallagher (2001) performed circular 
aperture photometry on MGG-F in F439W images from 
WFPC 2 and find a projected half-light radius of 160 ± 20 
mas (2.8 ± 0.3 pc). This value does not reflect deconvo- 
lution, but rather a "simple approach" of subtracting 150 
mas from the observed value in quadrature to correct for 
the PSF. By comparison, our profile-fitting method finds a 
smaller projected half-light radius of 89± 11 mas (1.6±0.2 
pc) in both F160W and F222M. Sternberg (1998) notes 
that if the massive red supergiants are spatially segregated 
from intermediate mass stars, the SSCs will appear more 
compact in the near-IR than in the visible. The smaller 
size of MGG-F as measured in the NICMOS images may 
thus be evidence of mass segregation. 

5. DISCUSSION 

The combination of spectroscopic and photometric data 
enables us to examine the stellar composition of the clus- 
ters. In particular, we can use measured light-to-mass ra- 
tios to draw inferences on the form of the cluster IMFs. 
Using the stellar velocity dispersions and half-light radii 
determined above, we find the cluster masses in §5.1. Ex- 
tinction to the clusters, which plays a critical role in sub- 
sequent analysis, is estimated in §5.2, and cluster ages are 
estimated in §5.3. The light-to-mass ratios and IMF impli- 
cations are discussed in §5.4, leading to discussion of the 
fate of the clusters in §5.5. 

5.1. Kinematic Masses 

The mass of a star cluster may be determined by ap- 
plication of the virial theorem, which dictates that for a 
cluster in equilibrium, kinetic and potential energy are re- 
lated by V 2 = rjGM/R. Spitzer (1969) shows that when 
expressed in terms of the radius containing half of the mass 
of the cluster (r^), the constant of proportionality is well 
approximated by the value n = 0.4 for most mass dis- 
tribution models. Assuming that the light profile traces 
the mass distribution, the half-light radius is used as a 
proxy for the half-mass radius. For an isotropic velocity 
dispersion, the half-light radius is related to the observed 
half-light radius in projection (r^ p ) by r^p ~ fr^ (Spitzer 
1987), and the three-dimensional rms velocity V is related 
to the measured one-dimensional line-of-sight velocity dis- 
persion oy by V 2 = 3d 2 . Armed with determinations of 
the velocity dispersion and projected half-light radius, we 
are thus able to estimate the mass of a virialized super 
star cluster: 

M = 10^^ 

Under these assumptions, we find a mass of 1.5±0.3x 10 6 
M© for MGG-9. This corresponds to an average mass den- 
sity of 1.1 ±0.3 x 10 4 Mq pc~ 3 within the half-mass radius. 
MGG-11 is less massive at 3.5±0.7x 10 5 M , with a higher 
average mass density within r^j, of 2.7±0.9x 10 4 M Q pc~ 3 . 

As noted in Section 4, MGG-F shows evidence of mass 
segregation. While extensive dynamical mass segregation 



in MGG-9 and MGG-11 is unlikely (Section 5.3), it is pos- 
sible that the most massive stars formed nearer the cluster 
cores. In this case, the red supergiant velocity dispersions 
we measure in the near-IR would be smaller than the clus- 
ter mean and the masses we derive would represent lower 
limits. 

5.2. Extinction 

As noted in Section 4, near-IR images of the nuclear 
starburst show substantial diffuse emission which hints at 
the importance of interstellar extinction (see also Figure 
6). Theoretical studies (e.g., Hills 1980) indicate that the 
star formation efficiency is high (> 50%) in massive star 
clusters. Considering the additional gas-dispersing effects 
of supernovae and stellar winds from young, high-mass 
stars, it is unlikely that much gas from the natal cloud 
remains in the clusters after a few million years. We have 
thus assumed a foreground screen model for our extinction 
estimates. 

The spectra of both MGG-9 and MGG-11 are dominated 
by the light of early M-type supergiants. Based on the in- 
trinsic color of an M2I star (representing the midpoint of 
the M0I - M4.5I range), we can estimate extinction along 
the line of sight to the clusters. We calculated synthetic 
magnitudes for stars in the Pickles (1998) spectral library 
for NICMOS filters. The Pickles library includes an M2I 
star, which has a synthetic [F160W] - [F222M] color of 
0.52. We adopt this as the intrinsic color of the clusters, 
and allow for a scatter of ±0.1 when calculating uncertain- 
ties. 

Evolutionary tracks from population synthesis models 
predict a similar color ([F160W] - [F222M] - 0.5) for the 
clusters, as seen in the color-magnitude diagram (CMD) 
for the nuclear starburst (Figure 7). The Starburst99 (Lei- 
thcrcr ct al. 1999) models used generate H and K mag- 
nitudes on the Johnson system. We used the method of 
Marleau et al. (2000) to determine the color transforma- 
tion between the Starburst99 output and the NICMOS 
filters: 

([F160W] - [F222M]) = 1.46(i? - K)j + 0.042 

Using the data in Table 3, the [F160W] - [F222M] col- 
ors of MGG-9 and MGG-11 are 1.35±0.16 and 1.07±0.16, 
respectively. Subtracting the adopted intrinsic color of the 
clusters, we find color excesses of E([F160W] - [F222M]) 
= 0.83 ± 0.19 for MGG-9 and 0.55 ± 0.19 for MGG- 
11. Applying the extinction law of Cardelli et al. (1989), 
we estimate the extinction in the F160W bandpass as 
A FW ow = 2.47£([F160W] - [F222M]). Thus we find 
Apieow = 2.1±0.5 for MGG-9 and 1.4±0.5 for MGG-11. 

As noted by Cardelli et al., there exists a single extinc- 
tion law for A > 0.90/im. As such, the estimated extinction 
Afwow is independent of assumptions regarding the value 
of Ry. For reference, if Ry = 3.1, ^4fi60W = 0.19^4y and 
if Ry = 5, Ai?i 6 ow = 0.22Ay. The clusters plotted in Fig- 
ure 7 display differential extinction, from zero to as high as 
perhaps 12 magnitudes in Ay. This is not surprising given 
the range of background intensities of the diffuse emission 
seen in the NICMOS images (Figure 1). Extinction plays 
a critical role in determining the light-to-mass ratios of 
the clusters (Section 5.4). Observers should be wary when 
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attempting to trace star formation with UBV photome- 
try and optical spectra in dusty starburst and interacting 
galaxies. 

5.3. Ages 

The ages of the SSCs may be constrained by several in- 
ferences drawn from their spectra. The absence of nebular 
emission lines indicates that there are no O stars remain- 
ing in the clusters. This sets a lower age limit of about 
6-7 Myr, the main sequence lifetime of O stars. As noted 
above, the spectra are dominated by the light of red super- 
giant stars as evidenced by the strong photospheric CO ab- 
sorption. This is broadly consistent with somewhat older 
starburst populations (ages ~ 10 7 yrs). 

Cool, intrinsically luminous asymptotic giant branch 
(AGB) stars contribute substantially to the near-IR light 
of clusters between ages 0.3 to 1.5 Gyr, dominating the 
near-IR light at ~ 600 Myr (Mouhcine & Langon 2002). 
The thermally-pulsing AGB (TP-AGB) phase begins to 
significantly affect the near-IR colors of a population at 
an age of about 100 Myr. The TP-AGB phase is not 
modelled by the Geneva evolutionary tracks used in Star- 
burst99, compromising the applicability of the population 
synthesis models beyond this point. The deep CO bands 
of TP-AGB stars could be mistaken for features of RSG 
stars in a younger population in the absence of other age 
indications (Gilbert 2002). Unfortunately, the best test 
for the predominance of AGB stars, the H2O/C2 index 
(Langon et al. 1999), falls just outside our spectral range 
at 1.77/xm. However, Forster Schreiber (1998) notes that 
none of the features expected of AGB stars are present in 
her near-IR spectra of the nuclear region of M82, setting 
an upper age limit of ~ 10 s yrs. 

We can estimate the age range more precisely through 
population synthesis modelling. Gilbert (2002) used Star- 
burst99 models to predict the flux-weighted average spec- 
tral type as a function of cluster age. Applying this tech- 
nique at _ff-band, we find that the light of a cluster popu- 
lation is dominated by M0 stars, as indicated by our cross- 
correlation analysis for these SSCs, at two different ages in 
a Salpeter IMF model: 5-7 Myr and 12-14 Myr. The for- 
mer can be ruled out in this there are no O stars 
present. During the intervening period (7-12 Myr), the 
dominant spectral type in f/-band is later than M0, but 
is at no time later than M3. The latest template star in 
our high-resolution NIRSPEC spectral atlas is type M4.5; 
this spectral type matches the cluster spectra as well as 
the M0 template. This suggests that the best match is 
likely M2 or M3, spectral types as yet missing from our 
template atlas. We therefore estimate ages of 7-12 Myr. 
This is consistent with the 10-12 Myr age of the "older" 
burst population found by Forster Schreiber (1998) for the 
nuclear starburst region. 

Figure 7 shows the CMD for the nuclear clusters. For 
context, we have ovcrplotted an evolutionary track based 
upon a model 10 6 M Q cluster with a Kroupa IMF (Section 
5.4) for stellar masses of 0.1 to 100 M Q . The clusters fall in 
a region which is consistent with ages of 8 to 60 Myr, and 
exhibit differential extinction (Section 5.2). Ages cannot 
be read directly from Figure 7 because of the dependence 
of F222M luminosity on total mass. Dereddening MGG-9 
until it meets the evolutionary track for a cluster with its 



derived mass of 1.5 x 10 6 M Q , we find a most likely age of 
approximately 10 Myr. When dereddened, MGG-11 is too 
bright for the Kroupa IMF evolutionary track correspond- 
ing to its derived mass at any age. Using the evolutionary 
track of a Salpeter IMF truncated at a minimum mass 
of 1 M Q , we find a most likely age for MGG-11 of ap- 
proximately 9 Myr. This is discussed further in Section 
5.4, along with predictions of the time evolution of light- 
to-mass ratios. It should be noted, however, that these 
age estimates are based on an assumed IMF, and cannot 
therefore be used to infer information on the stellar con- 
tent. For this purpose, we will use the IMF independent 
age estimates of 7-12 Myr. 

Measurement of the half-light radii and velocity dis- 
persions provides an estimate of the dynamical (crossing) 
times for the clusters. Both MGG-9 and MGG-11 have 
t cr <~ 10 5 years. The clusters are clearly hundreds of cross- 
ing times old and thus gravitationally bound. Calculating 
the characteristic relaxation time t r using the prescription 
of Meylan (1987), we find t r > 10 8 years for these SSCs. 
At the present young ages of the clusters, we do not expect 
significant, cluster-wide dynamical mass segregation. (See 
Section 5.4 for further discussion on mass segregation.) 

Smith & Gallagher (2001) estimated the age of MGG-F 
as 60±20 Myr. Our population synthesis models indicate a 
younger age of ~ 40 Myr. In the age range of 25-40 Myr, 
the cluster light in i?-band is dominated by the light of 
G-type stars. This implies E([F160W]-[F222M]) - 0.1 for 
MGG-F and thus Api60W ~ -2, in close agreement with 
the estimated E(B - V) = 0.9 ± 0.1 of Smith & Gallagher. 

5.4. Light-to-Mass Ratios and the IMF 

The near-infrared light of SSCs older than ~ 7 Myr 
is dominated by high-luminosity, massive, evolved stars. 
Lower-mass stars contribute very little to the spectrum 
of a young SSC, even if they comprise a majority of the 
cluster's mass. Measurement of the cluster's kinematic 
mass enables detection of the presence of low-mass stars 
independent of any assumptions regarding the form of the 
IMF. Quantification of the low-mass stellar content in turn 
provides insight as to whether the SSC can remain gravi- 
tationally bound over a globular cluster lifetime. 

From our photometry, we determine the total luminosity 
of the clusters in the F160W filter: L F160 = AncfF^Au 
(where Av = 1.37 x 10 13 Hz for the filter bandwidth of 
AA = 1177. 3A). Before adjustment for interstellar ex- 
tinction, we find Lfi 6 o/M of 0.28 ± 0.06 (in solar units, 
L /M Q ) for MGG-9 and 1.0 ± 0.2 for MGG-11. Dered- 
dening by the estimated extinctions (Section 5.2) gives 
L F i6o/M of 1.8^8 for MGG-9 and 3.5±?;| for MGG-11. 

Figure 8 shows predicted Lfiqq/M ratios func- 
tion of cluster age, based on population synthesis models. 
We have plotted models for three different IMF forms: a 
Salpeter (1955) power-law IMF with stellar masses of 0.1- 
100 Mq, a Kroupa (2001) IMF, and a Salpeter power law 
truncated with a minimum stellar mass of 1M Q . 

The stellar IMF is more complex than a single power 
law IMF extending down to the lowest stellar masses (see 
review in Meyer et al. 2000). In general, the IMF is repre- 
sented as a broken power law, with one slope for the distri- 
bution of high-mass stars and a lower slope or slopes for the 
low-mass (< 1M ) stars. The overall effect of such a dis- 
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tribution is to reduce the relative proportion of low-mass 
stars with respect to a single-slope power law distribution. 
The Kroupa IMF is respresentative: £(M) oc M~ ai where 

ai = 1.3, 0.1 < M/M© < 0.5 
a 2 = 2.3, 0.5 < M/M© < 100 

We find that the Lpi eo /M ratio of MGG-9 is consistent 
with a solar-metallicity Kroupa IMF for the cluster's age 
and dereddcned luminosity. It does not appear that there 
is any dimunition in the relative proportion of low-mass 
stars in this cluster. MGG-11, however, has a dereddened 
luminosity which is too bright for a Kroupa IMF at any 
age. MGG-11 appears to be deficient in low-mass stars. 
This suggests either a non-standard IMF or mass segrega- 
tion within the cluster (sec below). For comparison, Fig- 
ure 8 shows the Lpmo/M ratio of a single-power law IMF 
(with a — 2.35, as in a Salpeter IMF) which is truncated 
at a minimum mass of 1M©. Such a distribution of stellar 
masses would lead to a Lpi 60 /M ratio consistent with our 
observations for MGG-11. 

We also consider MGG-F. Using the mass estimate of 
Smith & Gallagher (2001) in combination with our pho- 
tometry, we find a dereddened Lpi^o/M = 0.4 ±0.1. This 
light-to-mass ratio is likely too high for the cluster's age 
and a standard Kroupa IMF, again suggesting a deficiency 
of low-mass stars. Smith & Gallagher conclude that the 
lower mass limit is 2-3M for a power-law IMF of slope 
a = 2.3; our lower estimate of the age suggests that the 
lower mass cutoff need only occur around 1M . 

As noted in Section 4, there is evidence of mass segrega- 
tion in MGG-F. The relaxation timescale is too great for 
significant dynamical mass segregation to have occurred 
throughout clusters MGG-9 and MGG-11 (ages <~ 10 7 yr). 
Numerical simulations by Bonnell & Davies (1998) indi- 
cate that although clusters younger than their relaxation 
timescale are not fully mass-segregated, they are likely to 
show some measure of mass segregation at young ages. 
Specifically, mass segregation of the most massive stars is 
known to occur more rapidly than relaxation of the entire 
cluster (Gerhard 2000). The ncar-IR light is dominated 
by RSG stars which are evolved from high-mass main se- 
quence stars — stars which are expected to approach en- 
ergy equipartition most rapidly. In addition, there is ev- 
idence of primordial mass segregation in young massive 
clusters such as the Orion Nebula Cluster (Hillenbrand & 
Hartmann 1998), R136 in 30 Doradus (e.g., Brandl et al. 
1996), and several clusters in the Large Magellanic Cloud 
(de Grijs et al. 2002b). If star formation processes have 
conspired to concentrate high mass stars at the center of 
the clusters, the relaxation timescale there will be shorter; 
dynamical mass segregation will proceed more rapidly at 
the core of the clusters, on the order of a few crossing 
times (de Grijs et al. 2002a). As a result, the cores of 
SSCs may undergo significant dynamical evolution in as 
little as 25 Myr (de Grijs et al. 2002c), somewhat younger 
than MGG-F. 

It is therefore not unlikely that there is significant mass 
segregation, either primordial or dynamical, in the clus- 
ters. In that case, the RSG starlight observed in F160W 

4 However, stars in the low velocity tail of the distribution function 
(e.g., Adams 2000). 



samples only the inner portion of the clusters which is 
dominated by evolved high-mass stars. The velocity dis- 
persions determined by our cross-correlation analysis are 
affected by only the mass located within the volume oc- 
cupied by RSG stars. If a significant amount of mass lies 
at larger radii, the half-light radius measured for the RSG 
stars will not trace the half-mass radius, leading to an un- 
derestimation of the cluster mass. 

If, for example, the half-mass radius of cluster MGG-11 
is greater than the measured RSG half-light radius by a 
factor of 1.5, the Lpieo/M ratio would be in line with a 
standard Kroupa IMF. Under such a segregation of masses, 
the RSG stars would occupy a volume filling factor of 30% 
within the cluster. 

5.5. Fate of the Clusters 

Whether SSCs MGG-9 and MGG-11 will resemble glob- 
ular clusters 10-15 billion years from now depends in 
large part upon two issues: whether the clusters can re- 
main bound despite mass loss due to stellar evolution and 
whether that evolution leads to a mass function simliar 
to old globular clusters. We will deal with these issues in 
turn. 

With ages in excess of 7 Myr, the OB stars formed in 
the clusters have begun evolving off the main sequence. 
We see no spectral evidence of nebular emission; stellar 
winds and ionization from the massive stars have cleared 
residual gas from the clusters' formation. In addition, any 
O stars in the clusters have likely exploded as supernovae, 
clearing any gas not removed gradually during the main 
sequence lifetime of the OB stars. Studies (e.g., Hills 1980; 
Goodwin 1997) have shown that loss of more than half of 
a cluster's mass in less than a crossing time can dissoci- 
ate the cluster. 4 The context for such studies is expulsion 
of residual gas and the implied star formation efficiencies 
of bound clusters. Although the residual gas has been 
removed from MGG-9 and MGG-11, the clusters will con- 
tinue to lose mass over time due to evolution of the mem- 
ber stars. If the mass loss is gradual, the stellar orbits 
have time to adjust to the new potential and the cluster 
will simply expand to a new equilibrium radius. The re- 
maining opportunity for substantial, impulsive gas loss is 
from supernovae of B stars. For a standard Kroupa IMF, 
14% of the cluster's mass is in stars with mass > 8M Q . 
Truncating the lower end of the IMF at 1M© shifts more 
of the cluster's mass to these supernova progenitors: 35% 
would be in stars with mass > 8M . In either case, how- 
ever, there is not enough mass for impulsive mass loss via 
supernovae to dissociate the clusters. 

While the clusters are likely to remain bound in spite 
of internal threats, they still must face various external 
destruction mechanisms (e.g., evaporation, disk-shocking 
and dynamical friction). Although the prospects of these 
clusters against these mechanisms is beyond the scope of 
the present work, we can make some general comments. 
Low-mass clusters are more easily dispersed by each of 
these processes (Ashman & Zepf 1998); the large masses 
of MGG-9 and MGG-11 (relative to Galactic globular clus- 
ters) will make these clusters less prone to disruption. The 
low galactocentric radii of the clusters, however, will tend 

remain bound in spite of loss of as much as 80% of the cluster mass 
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to make them more prone to disruption by dynamical fric- 
tion. 

The second issue is whether the SSCs will resemble glob- 
ular clusters after 10-15 Gyr of stellar evolution. The most 
metal-poor globular clusters in the Galaxy are about 15 
Gyr old (Vandenberg et al. 1996), with little or no dis- 
persion in age among the Galactic globular cluster system 
(Stetson et al. 1996). These old globular clusters have 
main sequence turn-offs at ~ 0.8 M , with higher-mass 
stars having evolved, leaving behind remnants. Sternberg 
(1998) quotes typical mass fractions for the component ob- 
jects in globular clusters: /(MS) - 0.73, /(WD) - 0.25, 
and /(NS) ~ 0.02 for the fraction of the mass contained 
in main sequence stars, white dwarfs and neutron stars, 
respectively. Applying Sternberg's simple model for map- 
ping main sequence progenitors to remnants, we are able 
to estimate the composition of a cluster after 15 Gyr of 
evolution (assuming the cluster remains bound and loses 
no stars, e.g. due to evaporation). 

For the case of a power law IMF truncated at a lower 
mass of 1 M Q , it is clear that no stars will remain on 
the main sequence after 15 Gyr; such a cluster could not 
resemble a globular cluster as described here. More inter- 
esting are the effects of evolution on the Kroupa IMF of 
Section 5.4. After 15 Gyr, the cluster will have lost 51% of 
its SSC mass (i.e., the mass present at an age of 10 Myr). 
The composition of the cluster would be /(MS) = 0.60, 
/(WD) = 0.37, and /(NS) = 0.04. This is similar to the 
typical globular cluster composition quoted by Sternberg, 
but with a larger proportion of the mass in the form of 
stellar remnants. Evolving MGG-9 for a globular cluster 
lifetime of 15 Gyr leaves a cluster of M= 8 x 10 5 M with a 
composition somewhat high in remnants but not dissimilar 
to Galactic globular clusters. 

Finally, we note that low-mass stars are preferentially 
lost in the process of stellar evaporation. This effect is 
exacerbated by mass segregation, as low-mass stars are 
in a shallower potential and have lower escape velocities. 
Dynamical effects are likely to decrease the mass fraction 
of low-mass main sequence stars, particularly in low-mass 
clusters for which these effects are strongest. 

6. CONCLUSIONS 

We have used high resolution near-IR spectra to mea- 
sure the velocity dispersions of two spatially-resolved su- 
per star clusters in the nuclear starburst of M82. We have 
combined these data with high-resolution NICMOS im- 
ages to determine the kinematic masses and light-to-mass 
ratios of these SSCs. The near-IR spectra of the clusters 
are dominated by the light of RSG stars, consistent with 
ages of 7-12 Myr. 

The derived masses, sizes and ages are consistent with 
the hypothesis that the SSCs are young globular clus- 
ters. Whether they can evolve into objects resembling old 



Galactic globular clusters depends upon the IMF of the 
clusters. When compared to population synthesis models, 
the light-to-mass ratio of MGG-9 appears to be consis- 
tent with a standard Kroupa IMF, and the cluster could 
resemble a globular cluster 15 Gyr from now. The light-to- 
mass ratios of MGG-11 and MGG-F suggest a deficiency 
of low-mass stars; it is unlikely that the composition of 
these clusters will resemble globular clusters after 15 Gyr 
of stellar evolution. There is evidence for mass segregation 
in MGG-F; if the red supergiant stars are sampling only 
the inner portion of the clusters, our mass estimates would 
represent lower limits. This could account for the apparent 
deficiency of low- mass stars and bring the Lpi eo /M ratios 
in line with standard IMFs. The presence of mass segre- 
gation would further suggest that the SSCs are scaled-up, 
high-mass analogues of young Galactic clusters (e.g., the 
Orion Nebula Cluster and NGC 3603). 

These results are critically dependent on extinction es- 
timates. The population of clusters in the nuclear star- 
burst suffers from a wide range of differential extinction, 
complicating optical studies. The lower extinction in the 
near-infrared is a distinct advantage in determination of 
cluster properties in young, dusty starburst regions. 

In order to assess whether the young SSCs in M82 can 
evolve into a population of old globular clusters, they 
should be characterized both as a system and as individ- 
ual objects. The present work represents a pilot project in 
preparation for a more complete study of the SSCs in the 
nuclear starburst. We have undertaken an observing pro- 
gram at the Keck Observatory to use NIRSPEC to mea- 
sure velocity dispersions for additional SSCs in the nu- 
clear starburst. With these spectra, we will determine the 
masses and light-to-mass ratios for the population of SSCs 
and examine the luminosity and mass functions of the clus- 
ter system. The presence of system- wide deviations from a 
standard IMF would be more compelling than suggestive 
results for individual clusters. 
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Fig. 1. — Mosaic of NICMOS F160W images of the nucleus of M82. Clusters included in Table 3 are labelled with the secondary identifier 
for reference. One of the four nod positions of the 24" X 0."432 slit is shown. At the distance of M82, 1" = 17.5 pc; the scale bar at lower 
left is 100 pc long. 
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Wavelength (/xm) 



Fig. 2. — Near-infrared rest-frame spectra of echelle order 46 for the clusters and the MOIab template star. The spectra arc identically 
normalized, but offset vertically for clarity (for reference, zero points are —0.9, 0.2 and 1.2, from bottom to top). Many of the features visible 
in the supergiant spectrum are evident in the SSC spectra, although they appear "washed out" by the stellar velocity dispersion of the cluster 
(e.g., the CO bandhead at 1.6618/im). The sky emission spectrum is plotted at the bottom for reference. 
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Fig. 3. — Near-infrared rest-frame spectra of echelle order 49 for the clusters and the MOIab template star. The spectra are presented in 
the same manner as Figure 2. This order is dominated by strong metal absorption lines, including Ti I at 1.5548 p.m. The CO 3-0 bandhcad 
at 1.5582 fim is also prominent. The sky emission spectrum is plotted at the bottom for reference. Pixels replaced by the local median value 
(Section 3) may be identified by comparison of the cluster and sky spectra. 
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Velocity (km/s) Velocity (km/s) 

Fig. 4. — Cross-correlation results for the clusters and the template MOIab star, averaged over all observed echelle orders (Section 3). The 
peak of the template star auto-correlation function (ACF) has been scaled to the peak of the CCF to emphasize the difference in width. 
The CCF is wider than the ACF due to the convolution with the velocity distribution function of the cluster, (a) CCF for MGG-9 and the 
template MOIab star HR 2289. The ACF peak of the template is clearly narrower than the CCF peak. The solid line is the antisymmetric 
part of the CCF, showing a slight asymmetry, (b) CCF for MGG-11 and the same template MOIab star, HR 2289. MGG-11 has a smaller 
velocity dispersion than MGG-9, as demonstrated by the smaller difference in the widths of the ACF and CCF. The CCF for MGG-11 is 
highly symmetric. 
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Fig. 5. — (a) Radial profile plot for MGG-9. The sky-subtracted NICMOS data (+) were fitted in two dimensions, but are collapsed 
azimuthally for this plot. The ordinate is plotted in [F160W] magnitudes per square arcsec (/1160). The solid line shows the azimuthal average 
of the fit, which is a King model convolved with the NIC2 PSF. The lower plot shows the two-dimensional, point-to-point residuals (data 
minus fit) in NIC2 counts/sec. (b) Radial profile plot for MGG-11. 




2468 2468 
arcsec arcsec 

Fig. 6. — A comparison of the NICMOS F160W image with an archival WFPC2 F814W image of the region surrounding clusters MGG-9 
and MGG-11. Circles mark the position of these clusters in each image; while there is faint emission evident at the position of MGG-11 in 
the F814W image, MGG-9 is completely invisible due to heavy extinction. The NICMOS image has been resampled to match the 0.1 arcsec 
pixels of WFPC2. Certain clusters have been labelled as reference points between the images; all axes are labelled in arcsec. 
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Fig. 7. — Color-magnitude diagram of the SSCs included in Table 2. MGG-9 (filled circle) is the reddest cluster in the sample, and MGG-11 
(filled square) is also more red than the average; MGG-F (filled triangle) is the bluest cluster. The Starburst99 evolutionary track for a 10 6 
Mq cluster with a Kroupa IMF is displayed with cluster ages labelled in millions of years. Note that the magnitude of a cluster is dependent 
upon its mass; the age of a cluster cannot be directly determined from the plotted track without knowledge of its mass. 
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Fig. 8. — Evolution of the luminosity-to-mass ratio as a function of age based on Starburst99 modelling. Three IMF models arc shown 
for a fiducial cluster mass of 10 Mq: Salpctcr singlc-powcr-law IMFs with minimum stellar masses of 0.1 or 1.0 Mq, and a Kroupa IMF 
(see Section 5.4). MGG-9 appears consistent with a Kroupa IMF, whereas MGG-11 and MGG-F may be deficient in low-mass stars. The 
vertical error bars for MGG-9 and MGG-11 demonstrate the uncertainty on the L/M ratio; their horizontal location does not reflect any age 
information within the range of 7-12 Myr. The MGG-F data are based on the mass estimate of Smith & Gallagher (2001). 
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Table 1 
Template Star Observations 



Template Star Spectral Type Date Airmass Atmospheric Calibrator 



HR 


3459 


G2Iab 


2000 


Dec 


10 


1.12 


HR 2432 (BHb) 


HD 


186293 


KOI 


2001 


Jun 


17 


1.08 


HR 7767 (09V) 


HR 


2119 


K0II 


2000 


Dec 


10 


1.26 


HR 2432 (BHb) 


HR 


8248 


Kllab 


2001 


Jun 


17 


1.11 


HR 7767 (09V) 


HD 


38054 


K3III 


2002 


Feb 


23 


1.53 


HD 32249 (B3V) 


HD 


201251 


K4Iab 


2001 


Jun 


17 


1.14 


HR 7767 (09V) 


HD 


185622 


K4Ib 


2001 


Jun 


17 


1.08 


HR 7767 (09V) 


HD 


200905 


K4.5Ib 


2001 


Jun 


17 


1.10 


HR 7767 (09V) 


HD 


43335 


K5Iab 


2002 


Feb 


23 


1.18 


HD 41753 (B3V) 


HR 


2289 


MOIab 


2000 


Dec 


10 


1.28 


HR 2432 (BHb) 


HD 


14488 


M4.5Iab 


2003 


Jan 


19 


1.26 


HD 14372 (B5V) 



Table 2 
Cross-Correlation Results 







MGG-9 




MGG-11 




Template 


Spectral 


Cross-Correlation 


OCCF° 


Cross- Correlation 


<JccF h 


Star 


Type 


Pcak a 


(km s- 1 ) 


Pcak a 


(km s _1 ) 


HR 3459 


G2Iab 


0.33 ±0.02 


17 ± 2 


0.33 ±0.03 


11±2 


HD 186293 


KOI 


0.42 ±0.02 


19 ± 1.2 


0.47 ±0.02 


11±2 


HR 2119 


K0II 


0.39 ±0.02 


20 ± 1.1 


0.44 ±0.02 


12 ± 1.1 


HR 8248 


Kllab 


0.62 ±0.01 


17.2 ±0.8 


0.69 ±0.02 


12.0 ±0.8 


HD 38054 


K3III 


0.52 ±0.01 


18.2 ±0.8 


0.59 ±0.02 


12.3 ±0.8 


HD 201251 


K4Iab 


0.59 ±0.01 


17.3 ±0.8 


0.66 ±0.02 


12.0 ±0.8 


HD 185622 


K4Ib 


0.68 ±0.01 


16.8 ±0.8 


0.75 ±0.01 


11.9 ±0.8 


HD 200905 


K4.5Ib 


0.64 ±0.01 


17.0 ±0.8 


0.72 ±0.01 


12.1 ±0.8 


HD 43335 


K5Iab 


0.58 ±0.01 


17.5 ±0.7 


0.67 ±0.01 


12.4 ±0.8 


HR 2289 


MOIab 


0.72 ±0.01 


16.2 ±0.8 


0.79 ±0.01 


11.3 ±0.7 


HD 14488 


M4.5Iab 


0.72 ±0.01 


15.5 ±0.7 


0.79 ±0.01 


11.5 ±0.8 



a Maximum value of the cross-correlation function. 

b o"ccF is the stellar velocity dispersion as determined from the cross-correlation function. 
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Table 3 

Super Star Cluster Parameters 



SSC 


MGG R.A.t Dec.t 


[F160IT]* 


[F222M]* 


Thp t 




(sec) (arcsec) 


(mags) 


(mags) 


(mas) 



J0955468+694937 


L 


46 


,8 


37 


,1 


12.05 


± 


0. 


,03 


11.45 


± 





.03 


83 ± 


10 


J0955471+694940 


F 


47. 


,1 


40. 


.8 


12.96 


± 


0. 


,04 


12.65 


± 


0, 


.04 


89 ± 


11 


J0955483+694941 


b 


48 


3 


41. 


.0 


14.5 


± 


0. 


,1 


13.46 


± 


0, 


.06 


113 ± 


14 


J0955488+694942 


c 


48. 


.8 


42. 


.5 


16.1 


± 


0. 


.4 


14.7 


± 


0. 


.5 


73 ± 


9 


J0955489+694940 


d 


48. 


.9 


40. 


.5 


15.5 


± 


0. 


,2 


14.5 


± 


0. 


,2 


94 ± 


11 


J0955502+694945 


11 


50. 


.2 


45. 


.1 


13.10 


± 


0. 


.15 


12.03 


± 


0. 


.05 


66 ± 


8 


J0955505+694945 


9 


50. 


.5 


45. 


.1 


12.90 


± 


0. 


,15 


11.55 


± 


0. 


.05 


146 ± 


18 


J0955506+694942 


8 


50. 


.6 


42. 


.9 


14.3 


± 


0. 


.3 


13.13 


± 


0. 


,15 


91 ± 


11 


J0955510+694946 


7 


51. 


.0 


46. 


.6 


13.4 


± 


0. 


.2 


12.80 


± 


0. 


.15 


156 ± 


19 


J0955513+694946 


6 


51 


.3 


46. 


.8 


13.03 


± 


0. 


,15 


12.18 


± 


0. 


,05 


79 ± 


9 


J0955525+694943 


k 


52 


,5 


43. 


.2 


14.1 


± 


0. 


,3 


13.8 


± 


0. 


,3 


166 ± 


20 


J0955526+694943 


m 


52. 


.6 


43. 


,0 


14.8 


± 





.7 


14.0 


± 


0. 


.3 


79 ± 


9 


J0955528+694950 


q 


52 


,8 


50. 


.6 


14.3 


± 


0. 


.3 


13.7 


± 


0. 


.2 


113 ± 


14 


J0955529+694944 


3 


52. 


.9 


44. 


.4 


13.3 


± 


0. 


.2 


12.37 


± 


0. 


.08 


90 ± 


11 


J0955535+694949 


la 


53. 


.5 


49. 


.7 


12.90 


± 


0. 


.15 


12.43 


± 


0. 


.08 


119 ± 


14 


J0955536+694948 


lb 


53. 


.6 


48. 


.4 


13.5 


± 


0. 


.2 


12.8 


± 


0. 


,1 


108 ± 


13 


J0955537+694949 


lc 


53. 


.8 


49. 


.4 


13.5 


± 


0. 


,2 


12.7 


± 


0. 


,1 


87 ± 


11 


J0955548+694953 


r 


54. 


.8 


53 


.9 


13.43 


± 


0. 


,06 


12.67 


± 


0. 


.04 


101 ± 


12 


J0955560+694952 


s 


56. 


.0 


52. 


.3 


15.25 


± 


0. 


.15 


14.2 


± 


0. 


.1 


105 ± 


13 


J0955571+694954 


t 


57. 


.1 


54. 


.0 


14.93 


± 


0. 


.15 


14.06 


± 


0. 


,09 


100 ± 


12 



t Offsets from coordinates: R.A. = 09 55 00.0, Dec. = 69 40 00.0 (J2000.0). The quoted 
astrometry is ±0."6. 

* Total magnitudes (uncorrected for reddening) based upon integration of King model 
fits to the data. 

■''The projected half-light radius, rh p , is based on King model fits to the F160W image. 



Note. - As recommended by the IAU, we have assigned each SSC a principal, 
coordinate-based designation. In addition, we refer to each source by a shorter nickname, 
prefaced by an acronym composed of the authors' initials, MGG. 



